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Loci /BLAST:  Overview 

►  Based  on  Loci  multi-physics  code  development 

framework 

°  Treats  applications  as  relational  databases  of  /acf5  (irregular  data) 
and  /Y//e5  (procedures  that  operate  on  the  data) 

°  Provides  automatic  code  generation,  parallelization,  dynamic 
memory  management  and  cache  optimization 

°  Lowers  cost  of  developing  compiex  multi-physics  applications 

►  Loci/CHEM:  Parent  of  Loci /BLAST 

°  Highiy-parallel,  full-featured  Eulerian  CFD  code  for  chemically 
reacting  flow  simulations 

°  Used  by  NASA,  USAF,  Army  and  Gov’t  iabs  aiong  with  Boeing  and 
several  small  aerospace  companies 

°  Extensive  V&V  using  MMS 
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Loci  Paradigm 


►  Define  computational  kernels  as  rules 

°  Framework  assembles  rules  into  applications  at  runtime  and 
schedules  execution 
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Loci  Rules  Example 

►  Loci  rules  mimic  mathematical  notation 

°  Parallelization  occurs  automatically  as  Loci  operates  on  collections 
of  entities 


//  Numerically  solve  the  partial  differential  equation:  ^=//V-Vi/ 
//  Using  Forward  Euler  time  integration 

//  Set  Initial  Value: 

$rule  pointwise(u{n=0}<-u_initial)  {  $u{n=0}  =  $u_initial  ;  } 

//  Forward  Euler  time-step:  =  t/”  H- Af/?" 

Srule  pointwise(u{n+l}<-u{n} ,R{n} ,dt)  {  $u{n+l}  =  $u{n}+$dt*$R{n}  ;  } 

//  Define  residual  (diffusion  operator) :  B  =  u  div(grad{u)) 

$rule  pointwise(R<-nu,div(grad(u) ) )  {  $R  =  $nu*$div(grad(u) )  ;  } 
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Loci  Multi-Physics  Framework 

►  Advantages: 

°  Al  detects  bugs  introduced  by  logically  inconsistent  specification 
of  numerical  models 

°  Decouples  computations  from  their  mappings  onto  distributed 
data  structures 

°  Databases  of  rules  become  knowledge  bases  for  composing 
increasingly  complex  simulations 

°  Parallelization  effort  is  greatly  simplified  and  is  automatic 

°  Optimizes  low  level  calculations  for  higher  levels  of  cache  use 
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Loci/BLAST  Capabilities 

►  Cell-centered,  finite-volume  method  for  general 
polyhedral  elements 

°  Overset  meshes  with  automated  hole  cutting 

►  2"^-orderTVD  Runge-Kutta  time  integration 

►  Approximate  Riemann  Fluxes  (HLLE,  HLLC) 

°  Robust  mixture  model  for  multi-material  flows 

►  Multiple  Equations  of  State 

°  Perfect  Gas 

°  Novel  tabular  EOS  based  on  Bezier  surfaces 
°  JWL  EOS  for  explosive  materials 
°  Linear  Barytropic  EOS  for  solids 
°  Multi-phase  EOS  for  soils 
°  Modified  Tait  EOS  for  water 
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Loci/BLAST  Capabilities 

►  Coupled  Lagrangian  particle  model  for  particulate 
flows 

►  Prescribed  burn  and  ignition  and  growth  reactive 
burn  models  for  explosive  detonation 

°  Secondary  combustion  model  for  non-ideal  explosives 

►  One-way  and  two-way  coupling  with  LS-DYNA  for 
Fluid-Structure  Interaction  (FSI)  simulations 

°  Conservative  load  transfer  between  Loci/BLAST  and  LS-DYNA 

°  Robust  deformation  of  CFD  mesh  in  response  to  structural 
deflection 

°  Simple  socket  protocol  provides  communication  between  the  two 
applications 
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Loci/BLAST  Soil  Model 

►  Soil  modeled  as  mixture  of  solids,  water  and  air 

°  Solids  composed  of  organic  mater,  clay  particles,  and  sand 

•  Currently,  solids  have  same  elastic  modulus  as  quartz 

•  No  material  strength  in  current  model 
°  Water  modeled  using  Modified  Tait  EOS 
°  Air  modeled  as  a  perfect  gas 

°  Mixture  mass  fractions  derived  assuming  that  the  soil  pore 
volume  can  be  determined  from  dry  soil 

°  Current  model  works  best  for  dryer  materials 

►  Soil  model  validated  using  mine  impulse  pendulum 
results 
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Underbody  Blast  Applications 

►  Loci/BI_AST-LS-DYNA  two-way  coupling  validated 
using  DRDC  Plate  experiment 

°  DRDC  plate  deflections  compared  to  results  of  Williams  et 
al  (7*^  internationai  LS-DYNA  Users  Conference,  2002)  for 
TACOM  Impulse  Loading  modei  (Westline,  1972) 

►  TARDEC  Generic  Hull  geometry  used  to  verify 
coupling  for  realistic  configurations 

°  Qualitative  results  only 
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DRDC  Plate  Experiment 


►  Large  weight  used  to  restrict 
target 

►  6kg  C4  at  5  cm  DOB 

►  31.75mm  thick  AL5083-H1  31 
target  plate 

►  Soil  density  of  2300kg/m3 


►  Various  soil  compositions 
tested 

►  Tabular  EoS 

►  Johnson-Cook  material 
strength  inputs 

►  4-noded  Belytschko-Tsay 
shell  elements 

►  0.5  cm  surface  resolution 


182  68  cm 


Mtine 


DRDC  Experiment  (Williams,  2002) 
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DRDC  Plate  Response 


Initial 


t=l  ms 


t=2ms 


Williams  (2002) 


Loci /BLAST  -  LSDYNA 
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DRDC  Plate  Response 
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TARDEC  Generic 

►  Geometry  represents  a 
notional  Army  vehicle 

►  Test  conditions: 

°  Charge:  6kg  cylinder  of  C4 

•  STANAG  4569  Level  2  mine 
blast  threat 

»  2  in  DOB 

°  Soil  taken  to  be  dry  sand  (70%  quartz  by  volume  fraction) 

►  Two  way-coupled  Loci/BLAST  -  LS-DYNA  analysis 

°  Conformal  meshes 

°  192  Loci/BLAST  processors,  1  LS-DYNA  processor 
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TARDEC  Generic  Hull 


Three  different  near-body  spacings 
used  in  CFD  mesh: 

Coarse  (5.0  cm  -  1 2M  cells) 
Medium  (2.5  cm  -  29M  cells) 

Fine  (1.25  cm  -  57M  cells) 

Soil  extends  3  ft  below  ground  plane 
Simulation  time  =  Sms 


Mesh 


Soil  volume  fraction  t=5  ms  Blast  front  t=5  ms 
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Mesh  Sensitivity:  Hull  Pressure 


Log^o  surface 
pressure  -1  ms 


P 

7.395e+00  r-j 
6.448C+00  m 
5.501C+00  ■ 
4.554C+00  ■ 
3.608e+00  H 
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Frame  and  Hull  Displacements 


0  0.001  0.002  Tinie(s|  0.003  0.004  0.005 
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Frame  and  Hull  Stresses  t=l  .5  ms 


Initial  frame  configuration 


Frame  stress  and  deflection 


iS-DVNA  keyword  deck  by  LS-PrePost 
Tim«  -  o,ooioo»^ 

Contours  of  Effootivn  Stress  iv-m} 
max  IP.  value 
mln=3illll,  at  «lain# 


Hull  stress  and  deflection 
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Conclusions  and  Future  Plans 

►  Loci/BLAST  is  a  modern  high-fidelity  blast 
simulation  tool  with  unique  features  not  found  in 
other  codes 

►  Successfully  applied  the  code  to  two-way  coupled 
underbody  blast  simulations 

°  Also  blast  in  urban  environments 

►  Future  plans 

°  Continue  verification  and  validation  efforts 
°  Add  adaptive  mesh  refinement 

°  Build  a  native  structural  dynamics  module  using  Loci  framework 
that  will  tightly  couple  with  CFD  components 
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Availability 

►  Loci/BLAST  and  Loci/CHEM  are  licensed  as  open  source 
software  and  are  available  free  of  charge.  They  are 
distributed  using  gpg  encryption  to  approved  users 
under  U.S.  ITAR  restrictions. 

»  http://www.simcenter.msstate.edu/software.php 

►  The  Loci  framework  is  open  source  and  available  from: 

°  http://www.cse.msstate.edu/~luke/loci/ 

►  Contacts 

»  Ed  Luke 

luke@cse.msstate.edu 
»  David  Thompson 

dst@cavs.msstate.edu 
°  Richard  Weed 

rweed@cavs.msstate.edu 
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Other  Loci  Applications 

►  Loci/CHEM:  Chemically  reacting  compressible  flow  solver. 

°  Currently  in  production  use  by  NASA  for  the  simulation  of  rocket 
motors,  plumes,  and  vehicles 

►  Loci/DROPLET:  Eulerian  and  Lagrangian  multiphase  solvers 

►  Loci/STREAM:  pressure-based  solver 

°  Developed  by  Streamline  Numerics  and  University  of  Florida 

►  Loci/FemLib:  Finite-element  linear  elasticity  thermal  stress  solver 
°  Developed  by  the  Cornell  Fracture  Group 

►  Loci/Radiation  model:  CA-DOM  non-gray  radiation  modeling 
°  Developed  in  collaboration  with  CFDRC 

►  The  Loci/THRUST:  High-Order  Discontinuous  Galerkin  Navier-Stokes 
Solver 

►  Various  multidisciplinary  simulation  tool  created  by  composing  the 
above  solvers 
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